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Figure 1: Measured [||] (symbols) and modeling spectra 
p(ui) of absolute motion and p(ui, L)/2 spectra of relative 
motion for the 2 AM SLAC site ground motion model. 

We believe that the frequency dependence of the mea- 
sured phase velocity v(f) is explained by the geological 
structure of the SLAC site where, as is typical, the ground 
rigidity and the density increase with depth. The surface 
motion primarily consists of transverse waves whose 
phase velocity is given by v s « y / E/(2p) and which are 
localized within one wavelength of the surface. If one 
plots the quantity v 2 /X versus wavelength A, we see that 
this value is almost constant, varying from 3000m/s 2 at 
A = 100m to 2000m/s 2 at A = 1000m. This is consistent 
with a ground density at the SLAC site that ranges from 
1.6 • 10 3 within the upper 100 m to 2.5 • 10 3 kg/m 3 at a 
kilometer depth and a Young's modulus E which increases 
from 10 9 Pa at 100 m to 10 10 Pa at 1000 m. These results 
seem to be quite reasonable for the SLAC geology, and, as 
we will see below, they also agree with explanations of the 
observed slow motion. 



3 SLOW GROUND MOTION 

Based on the arguments above, the wavelength at frequen- 
cies below 0.1 Hz quickly becomes much larger than the 
accelerator and eventually exceed the earth's size. In this 
case, the motion has little effect on the accelerator and at 
some point the notion of waves is not really applicable. 
Causes other than the wave mechanism must be responsi- 
ble for producing relative misalignments that are important 
at low frequencies. Such sources include the variation of 
temperature in the tunnel, underground water flow, spatial 
variation of ground properties combined with some exter- 
nal driving force, etc. These causes can produce misalign- 
ments with rather short wavelength in spite of their low fre- 
quencies. 

The ATL model of diffusive ground motion [Q] is an at- 
tempt to describe all these complex effects with a simple 
rule which states that the variance of the relative misalign- 
ment AX 2 is proportional to a coefficient A, the time T 
and the separation L: AX 2 = ATL. In the spectral repre- 
sentation this rule can be written as P(u>, k) = A/ (cu 2 k 2 ). 
It has been shown [|o[ that this rule adequately describes 
available measured data in many cases, however, typically 
only spatial or temporal information, but not both, was 
taken for a particular data set. Measurements where good 
statistics were collected, both in time and space and in a 
relevant region of parameter space, are sparse and difficult 
to perform. Thus, detailed investigation of slow motion is 
an urgent issue for future studies. 

The diffusive component of the ground motion model 
considered is based on measurements of slow motion per- 
formed at SLAC. First, measurements performed in the 
FFTB tunnel using the stretched wire alignment system 
over a baselength of 30 m give the value of A w 3 • 
10~ 7 /im 2 /(m-s) on a time scale of hours [B]. Second, a 
48 hour measurement of the linac tunnel motion performed 
with the linac laser alignment system over a baselength of 
1500 m gave A w 2 • 10~ 6 ^m 2 /(m-s) [|]. Finally, re- 
cent measurements using a similar technique were made 
over a period of one month and show that A w 10~ 7 - 
2 ■ 10~ 6 /im 2 /(m-s) for a wide frequency band of 0.01- 
10 _6 Hz [H], In the latter case, the major source of the 
slow 1/uo 2 motion was identified to be the temporal vari- 
ations of atmospheric pressure coupled to spatial variations 
of ground properties [B]. The atmospheric pressure was 
also thought to be responsible for a slow variation of the 
parameter A. 

The clear correlation of atmospheric pressure variation 
with deformation of the linac tunnel, observed in [|J, can 
only be explained if one assumes some variation of the 
ground properties along the linac. This variation can be 
due to changes in the Young's modulus E, changes in 
the topology such that the normal angle to the surface 
changes by Aa, or changes in the characteristic depth h 
of the softer surface layers. A rough estimate of the tunnel 
deformation due to variation of atmospheric pressure AP 
can be expressed as 
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Figure 2: Displacement of some points of SLAC linac tun- 
nel from 1966 through 1983 versus time and the approxi- 
mation in Eq. (|J) with r = 30 and to = 2 years. 
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The observed deformation of the tunnel AY = 50/im cor- 
responding to AP = 1000 Pa is consistent with this es- 
timation if AE/E~0.5, Aa -0.5 or Ah/h~0.5 and if 
one assumes E/h — 10 7 Pa/m. The former assumption 
is consistent with the heterogeneous landscape and geol- 
ogy at SLAC while the latter appears to agree well with the 
properties of the ground determined in the previous SLAC 
correlation measurements, if one assumes that h ~ A. 

No direct conclusions can be drawn from the measure- 
ments Q| to determine the spatial behavior of the observed 
slow motion because the relative motion was only mea- 
sured for one separation distance. However, the topology 
of many natural surfaces (including landscapes) exhibits a 
1/fc 2 behavior of the power spectra jilj]. Thus, it seems 
reasonable to expect that temporal pressure variation can 
also be a driving term of the spatial ATL-like motion. Fur- 
thermore, the measured parameter A can be extended from 
1500 m to a shorter scale, without contradicting the very 
short baseline measurements [||] which produced a similar 
value of A. 

It is also worth noting that the contribution to the param- 
eter A driven by the atmosphere scales as l/E 2 or as 
and therefore strongly depends on geology. Thus, the pa- 
rameter A, at a site with a much higher v s , would not be 
dominated by atmospheric contributions, while a site with 
softer ground and a v s half that at SLAC, may have a pa- 
rameter A as high as 3 • 10~ 5 /im 2 /(m-s). 

Finally, very slow motion, observed on a year-to-year 
time scale at SLAC, LEP, and other places, appears to be 




Figure 3: 17 year motion of the SLAC linac tunnel [hj] 



Figure 5: Rms relative motion versus time for L 
for the 2 a.m. SLAC site ground motion model. 
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systematic in time, i.e. AX Z 



r 2 



For example, 



measurements of the SLAC linac tunnel between 1966 and 
1983 10 show roughly linear motion in time with rates up 
to lmm/year in a few locations along the linac. Subsequent 
measurements indicate that the rate of this motion has de- 
creased over time although the direction of motion is still 
similar as is illustrated in Fig. ||. In the case of SLAC, the 
motion may have been caused primarily by settling effects, 
while in LEP, the cause may more likely be something dif- 
ferent such as underground water jl2t | . 

The temporal dependence of earth settlement problems 
typically are approximated as: 



1- 
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exp(-2.36t/r) (2) 



(1 + 2VVt) / 

where the typical value of t is years. This type of solu- 
tion exhibits \ft motion at the beginning which then slows 
and exponentially approaches Az/ max . An example of such 
a dependence is compared with the motion observed at 
SLAC in Fig. |. One can see that the early SLAC sys- 
tematic motion can be also described reasonably well by 
a linear in time motion, though nowadays the rate of the 
motion should be already much lower. 

The spatial characteristics of this systematic motion also 
seem to follow the 1/k 2 (or AX 2 oc L) behavior. This 
is evident in the displacements of the SLAC linac [[!]] af- 
ter 17 years which is shown in Fig. 0. The correspond- 
ing spatial spectrum is shown in Fig. [l] and it follows 
1/k 2 in the range of A from 20-500m. Although there 
is deviation from the 1/k 2 behavior at long wavelengths 



where there is limited data, this spectrum can be charac- 
terized as P sys t(t,k) « A sys tt 2 /k 2 with the parameter 



A sys t ~ 4 ■ 10 /xm /(m-s ) for early SLAC. An estimate 
of the rms misalignment due to this systematic motion is 
then AX 2 — A syst T 2 L. One can see that the transition 
between diffusive and systematic motion would occur at 
^trans = which in our case, assuming the value 

A = 5 • 10~ 7 /im 2 /(m-s) for the diffusive component of 
the SLAC ground motion model, would happen at about 



T tTi 
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The SLAC ground motion model includes all of the fea- 
tures that we have described. The transition from the 'fast' 
to the 'slow' motion is handled in a manner described in 
Ref. [ph. The absolute spectrum p(ui) and the spectrum 
of relative motion p(ui, L) are shown in Fig. |l|. The sys- 
tematic motion is not seen in this figure as it corresponds 
to much lower frequencies. However, it is seen in Fig. |] 
where the rms AX is calculated for L = 30 m by di- 
rect modeling of the ground motion using harmonic sum- 
mation [ 1 5 1 . One can see that this curve can be divided 
into three regions: wave dominated (T < 10 s), ATL- 
dominated (10^T^Sl0 5 s) and systematic motion domi- 
nated (T > T trans ~ 10 5 s). 

This ground motion model is included in the PWK mod- 
ule of the final focus design and analysis code FFADA [ fl3| ] 
which can perform analytical evaluations using the model 
spectra. The model is also included in the linac simulation 
code LIAR [ 14 1 where the summation of harmonics is used 



for direct simulations of the ground motion. 

4 CONCLUSION 

We have presented a model of ground motion for the SLAC 
site. This model includes fast, diffusive and systematic mo- 
tion with parameters that are consistent with the known ge- 
ological structure of the SLAC site. It is being now used to 
study the performance of the various systems in the Next 
Linear Collider. 

We would like to thank CAdolphsen, G.Bowden, 
M.Mayoud, R.Pitthan, R.Ruland, V.Shiltsev, and S.Takeda 
for various discussions of ground motion issues. 
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